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Abstract 

Pade approximants (PA's) are constructed from the perturbative coefficients 
of the free energy through 0{g^) in hot QCD. Pade summation is shown to 
reduce the renormahzation-scale dependence substantiahy even at tempera- 
ture (T) as low as 250 MeV. Also, PA's predict that the free energy does not 

deviate more than 10 % from the Stefan-Boltzmann limit for T > 250 MeV. 
PACS numbers: 12.38Cy, 12.38.Mh, ll.lO.Wx 
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Properties of quantum chromodynamics (QCD) at high temperature (T) acquire a lot 
of attention in relation to the physics of the early universe and of the relativistic heavy ion 
collisions [|1[. In non-abelian gauge theories, naive expectation, based on the asymptotic 
freedom, is that the perturbation theory works as far as the temperature is high enough [Q. 
However, this may not be the case because of the infrared sensitivity in the high-order terms 
beyond 0{g^) for the free energy with g being the QCD coupling constant [|]. 

Another practical problem is the asymptotic nature of the perturbation series, which 
can be seen from the expansion recently established through 0{g^) 0]. If one takes the 
expansion literally, the perturbation theory works only for T > 2 GeV which is order of 
magnitude larger than the critical temperature (Tc) of the QCD phase transition. On the 
other hand, the lattice QCD simulations for energy density, pressure and entropy density 
show that they approach the Stefan-Boltzmann (SB) limit rather quickly above — 150 
MeV [|. 

This poses a question that whether one could reconcile the non-perturbative lattice 
result with the perturbative expansion in the region < T < 2 GeV by making proper 
resummation of the perturbation series. From the experimental point of view, the highest 
temperature which one could access in the relativistic heavy-ion colliders are at most 300- 
500 MeV. Therefore, most of the future data reflect the physics in the above temperature 
interval. 

If the perturbative expansion has an asymptotic nature, the truncated series shows a large 
fluctuation for medium/strong couplings by the change of the order of the truncation. Never- 
theless, some information on the full result is reflected in the behavior of the fluctuation and 
the coeflicients of the perturbative expansion. This is the place where resummation method 
such as the Pade approach could play a role [0. In fact, the Pade summation has been suc- 
cessfully applied to the perturbative QCD series in high energy processes: it 'postdicts' the 
known higher-order terms, and also removes the renormalization-scheme dependence [pjQ. 

In the following, we will apply the Pade approach for free energy of QCD at flnite T with 
vanishing chemical potential and quark masses. The general structure of the free energy 
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F{T^ILi) for massless n/-flavors in MS scheme through 0{g^) reads 

Rif^) = = 1 + /2 fifi) + h fi^i) + (1) 

where /x is a renormahzation-scale, g = g/2iT and Fsb is the SB (classical) limit of the free 
energy. Since the explicit form of the coefficients fi is given in [|,|| , we will not recapitulate 
them here. 

We plot R{ii = 2tiT) forn/ = 4 in Fig. 1(a) as a function of ^^(/i = 2tiT) with as = 5'^/47r 
(see the discussions below for this choice of /x). The solid lines in Fig. 1(a) have a large 
fluctuation as one increases the order of the perturbation. as{27TT) used in Fig. 1(a) is shown 
in Fig. 1(b), where we take the two-loop /?- function with n/ = 4 for simplicity. Instead of 
expanding ag by l/ln(/x/A), we calculated the running of ag numerically with an initial 



condition ^^(/x = 5GeV) = 0.21 obtained in the MS scheme 



Before applying the Fade summation, let us first examine the renormalization-scheme 
dependence of i?(/x). Up to 0{g^) for the /^-function, the scheme dependence is equivalent 
to the renormalization-scale dependence [|n|]. The solid line with 'pert.' in Fig. 2 shows i?(/x) 
for n/ = 4 as a function of /x at T = 250 MeV. This temperature is rather low but is still 
1.7 times larger than Tc. In Fig. 2, i?(/i) shows a sizable /x dependence and even becomes 
unstable (negative pressure) for low /x. i?(/x) keeps this substantial /x-dependence unless T 
is extremely large. 

The principle of minimum sensitivity (PMS) criterion |ll[ does not work in the above 



situation, since there is no solution of the stability condition rfi?(/x)/rf/x = 0. The fastest 
convergence criterion (FAC) gives an unphysically large value /i = ST.TnT (for n/ = 4) 
The criterion motivated by the Brodsky-Lepage-Mackengie (BLM) [|1^] suggests /x = 
(0.95 ^ 4.4)7rT but gives rise to a series reliable only for T > 2 GeV Thus, any choice 
of /X cannot solve the problem at hand, and we really need a summation of the series. 

Perturbation series eq.(|T|) has two different features from that for high energy processes: 

(i) Odd powers of g appear. 



(ii) There arises In^' in the coefficients. In fact, the coefficient depends linearly on In^'. 



Since the standard Pade approximants are based on the ratio of polynomials, (ii) is a new 
feature beyond the standard method in a strict sense. In this paper, we take a simple 
procedure that In^ is regarded as a part of the coefficients 

Let us write down the general form of the Pade approximants (PA's) 



Cn and dm are the functions of Since we have 5 coefficients /i (1 < i < 5), we can 
construct PA's satisfying N + M <5. For example, the highest PA's are [3/2], [2/3], [4/1] 
and [1/4]. 

PA's discussed above are the most "naive" Pade approximants (NPA's). They have some 
unnatural feature from the physics point of view: NPA's contain terms proportional to g 
both in the numerator and the denominator. On the other hand, in eq.(|l]), we never have 
0{g) term because of a trivial reason (the free energy has no external legs). Although NPA's 
assure that their Tayler expansion by g does not contain 0(^)-term by the condition Ci = rfi, 
it is still unnatural that fictitious 0(^)-term play a crucial role to determine higher orders. 

Therefore, one may also try alternative PA's by assuming Ci = di = 0. In this "con- 
strained" Pade Approximants (CPA's), one can have PA's one-step beyond, i.e. N + M < 6. 
For example, the highest CPA's are [4/2], [2/4] and [3/3]. It turns out that the coeffi- 
cients (q, di) in CPA's are much simpler than those of NPA's because of the condition 
Ci = di = 0. In particular, (q, di) in NPA's usually become complicated functions of In^, 
while they simple depend linearly on In^ in CPA's. 

Besides NPA's and CPA's, one can also start with the effective charge aS(/x) = 
(i?(/x) — l)//2- The corresponding PA's for effective charge (EPA's) read aS'^/^1(/x) = 
f (1 + En=iCng'")/{l + E^=idmg'^). In this case, + M < 3 holds and highest PA's 
are [2/1] and [1/2]. However, one can show that [2/l]-EPA and [l/2]-EPA are equivalent to 
[4/l]-NPA and [3/2]-NPA respectively, so they are not independent from NPA. 



Also, it is expected to appear at 0{g^) level 0]. 



l + E: 




(2) 



4 



We have tried all possible PA's mentioned above and the following is a summary. 

(a) [3/2]-NPA, [2/3]-NPA, [4/2]-CPA, and [2/4]-CPA give results qualitatively consistent 
with each other. 

The sohd (dashed) lines in Fig.2 show the // dependence of CPA's (NPA's) at T = 250 
MeV. The renormalization-scale dependence is reduced substantially compared to the 
original series (the solid line with 'pert.'). This is one of the justifications that the Pade 
approximants point to the right direction. Similar /x independence has been reported 
on PA's for the Bjorken sum rules at high [§. In the latter case, the approximate 
invariance of PA's under the Euler transform is a key to the weak scale-dependence 
[|]. The similar argument is expected to hold at finite T. 

(b) [4/l]-NPA and [l/4]-NPA turn out to develop a pole in the denominator at ^ 0.04 
and 0.06, respectively. [3/3]-CPA also has poles at as ^ 0.12 and ^ 0.3. Whether 
they have real physical meaning or are the artifact of the approximation is not known. 

(c) Low-order PA's, namely + M < 4 for NPA's, + M < 5 for CPA's and A + M < 2 
for EPA's, do not properly 'postdict' the existence of In^ in /4 or the absence of In^ 
in /s. This implies that every known information of < 5) is required to construct 
proper PA's. 

Let us now focus on the successful cases in item (a) above. Because of the small 
dependence of PA's in (a), the choice of /x is not a serious problem anymore as far as 
/X > IGeV for T > 250MeV. The ELM scheme, which requires the leading Uf dependence 
to vanish, gives /x AAnT ( /x :^ ttT) when it is applied to the coefficient ( /s). Based 
on these observations together with the fact that 27rT is the lowest non-zero Matsubara 
frequency for gluons, we adopt /x = 2tiT as a typical renormalization scale. To see the 
quantitative prediction of the Pade summation, we plot, in Fig. 3, four different PA's for R 
as a function of ag with /x = 27rT. All four curves turn out to have similar trend: 
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First of all, they do not deviate more than 10 % from the SB limit even at T = 250 MeV. 
This is in contrast to 60 % deviation of R in eq.([l|) (see Fig. 1(a)). One should however note 
that such a small deviation does not necessary imply that the system can be approximated 
by non-interacting gas of quarks and gluons. There could be still strong interactions, but 
various effects tend to cancel for free energy. 

Secondly, PA's for the pressure (= — free energy) in Fig. 3 show that overall attraction 
(repulsion) occurs at high (low) T by the quark-gluon interactions. 

To compare PA's with the original series, we plot the result of [4/2]-CPA by the dashed 
line in Fig. 1(a) as an example. This shows the above two aspects explicitly. 

Our resummation is solely based on the perturbative coefficients calculated through 
0{g^) where the non-perturbative magnetic-mass of 0{g'^T) does not play a role. Therefore, 
prediction of PA's in this paper could be substantially modified in the strong coupling (low 
T) regime. Nevertheless, it will be interesting to estimate the 0{g^) coefficient by our PA's. 
The formulas become particularly simple for the cases [3/2]-CPA and [2/3]-CPA, in which 
/e = hh/fs and fe = /| + /| + {f^ - fl)/ h are obtained, respectively. In both cases, /e 
has a proper structure, namely a + 61n^, since only depends linearly on In^. This Pade 
prediction will work only for a part of /e which is insensitive to the physics of 0{g'^T). 

One should also note that we have not attempted to study the convergence property of 
PA's. This is because low-order PA's do not 'postdict' the known perturbative coefficients 
and thus do not capture the essential features of the full result as we mentioned. 

In summary, we have examined Pade approximants of the free energy of hot QCD. Naive 
PA's as well as constrained PA's reduce the renormalization-scale dependence substantially 
even at temperature as low as 250 MeV. After the Pade summation, the free energy takes 
a value close to the SB limit even at 250 MeV. Similar analyses with finite chemical poten- 
tial and finite quark masses, which are more relevant to the data in relativistic heavy-ion 
collisions, remain as future problems [1^ . 



After the submission of this work, the author was informed about a recent paper treating 
NPA's in QCD and 0^ theory (B. Kastening, |hep-ph/97082T9D . CPA's and EPA's are not 
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studied in this paper. 
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Fig.l: (a) Normalized free energy R{27rT) as a function of the running coupling as{27rT). 
0{g^) (/ = 2, 3, 4, 5) shows R summed up to the l-th order. The dashed line is the result 
of [4/2] Pade approximants. (b) as{27rT) as a function of T numerically calculated 
with the two-loop /^-function. 

Fig.2: Renormalization-scale dependence of the normalized free energy The solid line 

with 'pert.' corresponds to the perturbative evaluation up to 0{g^). Others are the 
results after Pade summation. 

Fig.3: Pade approximants of R{2tiT) as a function of as{2TiT). 
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